METAL-INSULATOR TRANSITION IN THE TWO-DIMENSIONAL HUBBARD 
MODEL AT HALF-FILLING WITH LIFETIME EFFECTS WITHIN THE 

MOMENT APPROACH. 

J.J. Rodriguez- Nunez 
Departamento de Fisica - CCNE, 
Universtdade Federal de Santa Maria, CAMOBI, 
97105 - 900 Santa Maria - RS, Brazil, 
e-m: jjrn@mfis.ccne.br 

S. Schafroth 
Physik-Institut der Universitat Zurich, 
. . . Wmterthurerstrasse 190, CH-8057 Ziirich, Switzerland 

' e-m; schafrot@physik.unizu.ch 

ON ! (February 1, 2008) 

We explore the effect of the imaginary part of the self-energy, ImT,{k,u)), having a single pole, 
n(fc,aj), with spectral weight, a(k), and quasi-particle lifetime, r(A:), on the density of states. We 
D ' solve the set of parameters, Q,{k,uj), a.{k), and T{k) by means of the moment approach (exact sum 

^0 ' rules) of Nolting Our choice for E(fc, w), satisfies the Kramers - Kronig relationship automati- 

cally. Due to our choice of the self - energy, the system is not a Fermi liquid for any value of the 
interaction, a result which is also true in the moment approach of Nolting without lifetime effects. 
By increasing the value of the local interaction, U /W , at half-filling (p = 1/2), we go from a para- 
^ . magnetic metal to a paramagnetic insulator, (Mott metal - insulator transition {MM IT)) for 

c 2 ^ ' values of U /W of the order of U /W > 1 {W is the band width) which is in agreement with numerical 

00 , results for finite lattices and for infinity dimensions (D — oo). These results settle down the main 

■ weakness of the spherical approximation of Nolting |l| : a finite gap for any finite value of the inter- 

action, i.e., an insulator for any finite value of U/W. Lifetime effects are absolutely indispensable. 
Our scheme works better than the one of improving the narrowing band factor, B{k), beyond the 
spherical approximation of Nolting [tl| . 
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After the discovery of the high- Tc materials js), the study of correlations has gained interested due to the fact 
that there is the behef Q that the normal properties of these materials could be explained in the framework of the 
Hubbard model |^,D, since electron correlations are strong, i.e., the on-site electron-electron repulsions U are much 
larger than the energies associated with the hybridization of atomic orbitals belonging to different atoms 0. The 
strategy, according to Anderson Q to address the problem of the high-Tc superconductivity, is to try to find a theory 
accounting for the normal state properties of the cuprates and then find an electron mechanism which destabilizes 
the normal state towards a superconducting state. The Hubbard Hamiltonian is a kind of minimum model |^ which 
takes into account quantum mechanical motion of electrons in a solid, and nonlinear repulsion between electrons. 
Even thought this model is too simple to describe solids faithfully, serious theoretical studies have revealed that to 
understand the various properties of it is a very difficult task, since is the simplest many-body Hamiltonian one can 
write down and which cannot be reduced to a single-particle theory px| ]. Its study will prove useful in developing 
various notions and techniques in statistical physics of many degrees of freedom. Besides cuprate superconductors, the 
Hubbard model may also be applicable to describing the metal - insulating transition in materials like LuxSri-xTiO^ 
and V2-J/O3 in which paramagnetic metal, antiferromagnetic insulator, and a phase antiferromagnetic metal are clearly 
separated in temperature - pressure phase diagram pT| |. All these features, found in the Hubbard model, make it 
suitable to reproduce experimental data. Additionally, the study of correlations by itself in the Hubbard model is a 
rewarding task since it will shed light on still unsolved points of the novel materials. For example, at high temperatures 
(Tc 30 — 130 K) these HTSC cuprates, which are poor conductors, become superconductors. This feature is strange 
indeed because the Coulomb repulsion is strong. Furthermore, the behavior of these materials at T > Tc is even 
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more puzzling than the superconductivity itself. Contrary to the predictions of the Fermi liquid theory, the resistivity 
at T > Tc and optimum doping is linear in temperature, i.e., R ^ T. This suggests a very strong scattering of 
elementary excitations. A discussion of the possible breakdown of Fermi liquid theory is given in Ref. fl^ . There is 
also another approach which consists in studying toy models, i.e., exactly solvable, [ p^ in order to get some ideas on 
the one-particle properties of highly correlated electron systems. 

In this Letter, we will use the moment approach (or sum rules) of Nolting for the spectral density, A{k,Lu). As 
it is well known in the literature fl^ j, one of the drawbacks of the moment approach in the spherical approximation 
- where B{k), the narrowing band factor, is not k - dependent - is that we always find a gap in the density of states 
{DOS). If the chemical potential happens to be in this gap, then we always have an insulator. It has been argued 
that the way to cure this unrealistic gap is to have a better approximation for the narrowing band factor, B{k), task 
which has proven to be hard just to the present 1 14|. It is due to that we have followed a different path which consists 
in proposing a single pole structure in the self-energy, J^{k,uj). 

The model we study is the Hubbard Hamiltonian 

H = t^^ci H n^- — act , (1) 

where cl (c^^.) are creation (annihilation) electron operators with spin a. n^^ = cl c^^. U is the local interaction, 
/d the chemical potential and we work in the grand canonical ensemble. We have adopted Einstein convention for 
repeated indices, i.e., for the Ng sites z, the z nearest-neighbor sites (n.n.) j and for spin up and down (cr = — = ±1). 
iij = —t, for n.n. and zero otherwise. 

Let us propose for the self - energy, I]{k,uj), the following single pole Ansatz 

o.A^^ -T^.A = "(fc),r(fc)G5R (2) 

LU — lt[k) — «I (k) 

With our choice for E(fc, w), we have the real and imaginary parts of the self - energy satisfy the Kramers - Kronig 
relations ||l5[| . 

By definition the one-particle Green function, G{k,uj), in terms of 'Z,{k,uj) is given as 

G(fc,w) = ^ , (3) 

LU - e^ - l^[k,uj) 

where = —2t{cos{kx) + cos{ky) — fi. Also, we will require the spectral density, A(fc, w) which is defined as 

A(k,u!) ^ -- lim ImG{k,uj + iS) (4) 

TT (5^0+ 

Using Eqs. (|| - ^), we arrive to the following expression for the spectral density 

A(fc»^zi (5) 



Using the first three sum rules of Nolting for the the spectral function of Eq. (|j) we obtain the following 
equations 

-t-oo 

A{k, uj)du; — ao{k) 

-oo 
4-00 

LuA(k,Lu)duj — ai(fc) 

-oo 

+ 00 

2 



uj A{k,uj)dLLj — a2{k) , (6) 

where the ai(fc)'s, i = 0,2 are given in Ref. Hi (see also Ref. ||). We do not use the forth moment or sum rule because 
we are at half - filling. Then, we assume that at p = 1/2 the chemical potential, /i = U/2. This is verified as we 
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will see from our results. At this moment, we would like to say once more that the drawback of previous calculations 
(included ours) consists in that the density of states which results of the two pole Ansatz for the one - particle Green 
function, in the spherical approximation of Nolting ffl, always has a gap. This solution (always a gap) is known in 



the literature as the Hubbard-I solution |l6l which has been critized since many years ago by Laura Roth |17| among 
others. We call the attention to Ref. |T8| where the authors point out to the fact that the fc-dependence has to be 
included in B{k). However, they study the negative Hubbard model in the strong coupling limit and in this limit we 
have a well developed correlation gap, anyway. Oles and coworkers are practicioners of the two - pole ansatz for the 
one - particle Green function, without lifetime effects. Here we are including lifetime effects as a crucial ingredient in 
the formulation. Morover, our results clearly show that the cure the Hubbard-I solution in the moment approach is 
not in the B{k) - term. 

In Fig. 1 we present the energy spectrum of the self - energy for several values of the Brillouin zone, for U /t — A. 

In Fig. 2 we show the fc-dependence of a{k) for U /t ~ 4, 6, 8, 12. In Fig. 3 we present the k dependence of r(fc) for 
several values of the Brillouin zone, for J7/t = 4, 6, 8, 12. In Fig. 4 we show the density of states, N{lo) vs lo. N{lo) is 
defined as 

* k 

In all the figures we have worked with a square lattice with periodicity of A^^ = 32 x 32. We conclude that the spectral 
weight, a{k) and the damping factor, r(fc), do not strongly depend on fc, for small values of U/W. For example, for 
U/t — 4.0, r{k)/t w 1.5. However, for larger values of U/W, a{k) has some fc dependence. On the contrary, r(fc) 
shows a stronger fc - dependence for larger values of U/W. We also note that a{k) x r(fc) < in order to give a 
positive spectral density (see Eq.(^)). Taking a look to the density of states, A^(a;) vs w, we see that the correlation 
gap opens up for U/W > 1.0. This is equivalent to the Hubbard-III like solution Due to our choice of the self 

- energy, we do not have a Fermi liquid . Edwards and Hertz po| have studied the breakdown of Fermi liquid 
theory in the Hubbard model at T = 0. They have drawn a phase boundary between a Fermi liquid and a non - 
Fermi liquid (p vs U in the paramagnetic phase). According to these authors for small values of U we have a Fermi 
liquid behavior (metallic Fermi liquid) and for larger values of U we have a non - Fermi liquid metal. In our approach 
we have a non - Fermi liquid behavior for any value of interaction. In order to agree with the authors of reference 
pO| we have to include a ui^ behavior for frequencies close by to the chemical potential in addition to the single pole 
structure chosen in the present work. A metallic Fermi to metallic non - Fermi liquid behavior has been obtained by 
Figueira, Anda and Nogueira pi]] . We should mention that the phase diagram of the Hubbard model contains the 
anti - ferromagnetic transition | |22| which has not been considered here. 

In short, we have postulated a one - pole Ansatz for the self - energy (Eq. (^)) in the moment approach of Nolting 
[01 . Our essential new working idea with respect to the now canonical method of Nolting is that we start from the 
self - energy while Nolting proposes a two pole ansatz for the one - particle Green function from the begining. In 
both approaches sum rules for the spectral function, A(fc,w), are imposed. As function of interaction, we see a Mott 
insulator - transition (MMIT) from a metal at weak interaction (no gap at the zero frequency) to a well developed 
gap (an insulator) for large values of the local interaction. Our formulation is not equivalent to the infinite dimension 
calculation of Georges et al |2^ since we have included the fc - dependence of the selfenergy (we are working in two 
dimensions). There is the calculation of Figueira, Anda and Nogueira |^ where the authors have also neglected the fc 

- dependence of I](fc, w). In particular, the authors of Refs. |^,^ find a Kondo peak at a; = 0. In our case, p = 1/2, 
the peak-like structure seen at the chemical potential is most likely due to the van Hove logarithmic singularity which 
always goes away with the opening of the correlation gap. With our precision p^ , we cannot conclude about the 
presence or not of the Kondo peak. Nolting himself has also studied the effect of damping on magnetism. We 
leave for the future the comparison with Ref. |2^ where the authors discuss the use of the moment approach for 
interpolating S(fc,a;) between weak and strong interaction. We include Figure 5 to present ri(m,n) vs n for fixed 
values of m, for U/t = 4.0 and a mesh of 10 x 10. It is worth looking to the band structure for m = 0. This curve has 
quasi the same structure than the band structure shown in Figure 1, for a bigger mesh (32 x 32). Then, we can say 
that finite size effects are minimal at least with respect to the band structure calculations. (See the discussion of Ref. 
p^). Finally, we would like to say that the non Fermi liquid behavior that we have in our approach (by construction) 
has nothing to do with any microscopic mechanism. Our arguments rely on general grounds of a self - energy ansatz 
which is used in the one - particle Green function to which exact sum rules (moments) are imposed. Those interested 
in physical mechanisms of non - Fermi liquid behavior are encouraged to take a look to the nice account presented in 
Anderson's book [^ . 
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FIGURES 



Figure 1. Oo,n vs n for U/t = 4.0 and p = 0.5. The k vector is given by fc = 7r(0, n)/32. As we work in two dimensions, 
the bandwidth isW = 8t. Compare with m = of Figure 5. 

Figure 2. a(0, n) vs n for different values oi U/t. Same parameters as in Figure 1. 

Figure 3. We show r(0,n) vs n for different values of U/t. Same parameters as in Figure 1. 

Figure 4. N{lo) vs to for U/t = A, 12. We see the opening of the correlation gap for U/W > 1.0, signaling the Mott 
metal - insulator transition. 

Figure 5. ^m,n vs n for some fixed values of m. U/t = 4.0, p = 0.5 and our mesh is 10 x 10. Compare m = with 
Figure 1. 
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